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Abstract 

An  analysis  has  been  developed  to  model  a  rotorcraft  with  power  loss.  Partial  power  recovery  to  a  flyaway  condition,  partial 
power  landing,  and  autorotation  to  landing  with  a  complete  power  loss  are  considered.  An  optimal  control  procedure  is 
implemented  to  predict  the  pilot  response  to  the  power  loss  with  the  objective  of  minimizing  the  horizontal  and  vertical 
impact  velocities  in  the  case  of  a  landing  and  achieving  a  sustainable  steady  state  flight  condition  for  flyaway.  The  physical 
and  optimal  control  models  are  based  on  previous  work,  but  have  been  extended  with  rate  controls  and  additional  constraints 
to  more  closely  represent  the  pilot  reaction  and  physics  of  the  helicopter  motion.  The  model  is  validated  for  autorotative 
landings  from  hover  and  forward  flight  initial  conditions  with  data  from  a  high  energy  rotor  system  flight  test.  With 
appropriate  inputs  to  simulate  the  flight  paths  in  the  test  data,  good  agreement  is  observed  between  the  optimal  control 
solutions  and  the  test  data.  An  example  flyaway  case  is  presented  to  illustrate  the  strengths  and  weaknesses  of  the  current 
implementation. 


Notation 

A  rotor  disk  area 

Cd0  rotor  blade  minimum  drag  coefficient 

Cj  rotor  thrust  coefficient 

Ct0  rotor  initial  thmst  coefficient  (weight) 

(Cj / cj )s  rotor  stall  thrust  coefficient 

Dx  vehicle  horizontal  drag 

Dz  vehicle  vertical  drag 

&  optimal  control  integral  part  of  objective  func¬ 

tion 

A  optimal  control  non-integral  part  of  objective 

function 

JA  optimal  control  objective  function 

Pk/d  Probability  of  kill  given  damage 

Q  rotor  torque 

R  rotor  radius 

S  optimal  control  non-differential  constraints 

T  rotor  thrust 

U  optimal  control  control  variable 

Ur )  optimal  control  slack  control  associated  with 

0 

Wx  optimal  control  velocity  weighting  factor 

X  optimal  control  state  variable 

a  rotor  blade  lift  curve  slope 

/  flat  plate  drag  area 

// ,  f(j  induced  velocity  and  ground  effect  factors 

go  acceleration  due  to  gravity,  non-dimension- 

alized  by  Q.qR2 
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ns  rotor  stall  exponent 

nt  number  of  time  steps  in  optimal  control  solu¬ 

tion 

Vcrit  vertical  impact  velocity  above  which  impact 

is  considered  an  attrition 
vx ,  V-  horizontal  and  vertical  velocities 

x,z  coordinate  system  horizontal  and  vertical  di¬ 

rections  and  displacements  (positive  forward 
and  down) 

fin  initial  rotor  speed 

II  optimal  control  scalar  parameter 

a  rotor  disk  angle 

y  rotor  Lock  number 

A  rotor  inflow 

p  rotor  advance  ratio 

co  non-dimensional  rotor  speed 

p  air  density 

a  rotor  solidity 

T  non-dimensional  time  (=  t£l o) 

T f  non-dimensional  terminal  time  (=  tfClo ) 

cp  optimal  control  differential  constraints 

t f/  optimal  control  end  constraints 

()v  differentiation  with  respect  to  normalized, 

non-dimensional  time 

()  *  differentiation  with  respect  to  non-dimension¬ 

al  time 

()o,()|o  initial  quantity 

Qt, ()|i  final  quantity 

( )  scale  factor,  specifically  for  X,  U,  and  II 

( )  slack  variable  associated  with  ( ) 

( )  min ,  ( )  max  minimum  and  maximum  values  for  a  variable 
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Fig.  1.  Example  Pk/d  plot  for  power  off  landing  showing 
areas  of  attrition  at  low  speed  and  forced  landing  at  high 
speed. 


Introduction 


The  previous  work  was  modified  and  extended  in  several 
ways  for  the  current  effort.  A  basic  departure  was  the  choice 
of  variables  for  the  equations  of  motion.  In  references  1-3, 
the  helicopter  was  controlled  by  thrust  coefficients  in  the 
horizontal  and  vertical  directions.  For  the  current  work,  a 
single  thrust  coefficient  and  disk  angle  were  used  instead. 
The  angle  introduced  nonlinear  sine  and  cosine  functions 
into  the  equations,  but  provided  several  advantages. 

In  the  previous  work,  only  a  complete  engine  failure  re¬ 
sulting  in  autorotation  was  addressed.  The  current  effort 
includes  partial  power  both  for  landings  and  mission  abort 
flyaways. 

Even  with  these  changes,  when  exercised  over  a  broad 
range  of  forward  speeds  and  initial  altitudes,  beyond  what 
was  considered  by  the  previous  researchers,  several  defi¬ 
ciencies  were  revealed  which  resulted  in  optimal  control 
paths  which  were  not  physically  possible.  These  are  ad¬ 
dressed  in  the  current  work,  mostly  by  additional  con¬ 
straints  on  the  flight  path,  described  later. 

In  this  paper,  a  description  of  the  initial  model  patterned 
after  previous  work  is  provided.  Following  that,  the  addi¬ 
tional  modifications,  constraints  and  the  justifications  for 
them  are  provided. 


The  purpose  of  the  present  work  is  to  develop  a  tool  to  be 
used  in  assessing  the  survivability  of  rotorcraft  with  power 
loss.  Specifically,  given  an  initial  altitude,  forward  speed, 
and  power  level,  the  tool  is  to  determine  if  the  helicopter 
can  fly  away  or  not,  and  if  not,  the  level  of  damage  that  will 
result  from  the  impact.  The  analysis  is  run  over  a  range 
of  altitude  and  forward  speed  and  the  results  are  collected 
into  a  probability  of  kill  given  damage  or  l\u.  The  Pk/d 
indicates  the  likelihood  of  each  type  of  kill  for  the  failure 
within  a  range  of  altitude  and  forward  speed. 

Three  kill  types  are  considered.  For  a  partial  power  loss, 
all  three  types  of  kill  are  possible.  First,  there  is  the  poten¬ 
tial  for  a  mission  abort,  where  the  helicopter  can  fly  away 
and  return  to  a  friendly  base  under  its  own  power.  If  the 
helicopter  is  unable  to  fly  away,  it  must  descend  to  landing 
and  potentially  sustain  additional  damage.  Only  two  levels 
of  damage  are  considered  for  landings,  namely  forced  land¬ 
ing,  where  the  helicopter  can  be  recovered,  repaired,  and 
returned  to  service,  and  attrition,  where  the  vehicle  is  a  to¬ 
tal  loss  and  is  removed  from  the  inventory.  Clearly,  for  a 
complete  power  loss  only  the  latter  two  kill  types  are  pos¬ 
sible.  An  example  Pk/d  plot  for  power-off  landing  is  shown 
in  Figure  1 . 

Further  objectives  for  the  project  were  that  it  be  com¬ 
putationally  simple  enough  to  run  on  a  workstation  com¬ 
puter  and  not  require  extensive  mathematical  or  engineering 
knowledge  of  the  user.  Given  those  objectives,  the  effort 
is  based  upon  and  extends  previous  work  from  the  1970’s 
and  1980’s,  where  similar  analysis  was  successfully  imple¬ 
mented  for  hover  (Ref.  1)  and  forward  flight  (Refs.  2, 3). 


Optimal  Control  Model 

The  bulk  of  the  analysis  is  an  optimal  control  procedure 
which  calculates  the  pilot  inputs  and  the  resulting  flight  path 
from  those  inputs.  The  algorithm  is  called  the  Sequential 
Gradient-Restoration  Algorithm  (Ref.  4).  In  the  algorithm, 
an  optimal  path  is  found  to  minimize  an  objective  function, 
consisting  of  an  integral  portion,  &  and  a  non-integral  por¬ 
tion  Sf.  The  solution  is  subject  to  differential  constraints  (p, 
non-differential  constraints  S,  and  end  constraints  if/  over  a 
unit  length  of  time.  The  constraints  are  specified  in  equa¬ 
tions  that  sum  to  zero  when  the  constraint  is  met.  If  the 
equation  is  nonzero,  the  value  is  constraint  error.  The  al¬ 
gorithm  operates  on  state  variables,  control  variables,  and 
scalar  variables  to  reduce  constraint  error  to  zero  and  to 
minimize  the  objective  function  as  much  as  possible.  It  is 
general  purpose  and  then  entire  problem  is  defined  by  the 
variables  described  above. 

Each  state  and  control  variable  is  a  vector  of  length  nt. 
The  scalar  parameters  are  values  which  are  variable  but  do 
not  change  with  time.  For  the  present  application,  the  only 
scalar  parameter  is  t,  which  is  the  scale  factor  between 
the  unit  optimization  interval  and  the  non-dimensionalized 
time.  Reference  4  provides  an  extremely  detailed  descrip¬ 
tion  of  the  algorithm  and  example  problems  with  tabulated 
data  to  verify  that  the  algorithm  is  implemented  correctly. 
Only  enough  explanation  is  provided  here  to  provide  con¬ 
text  for  the  application  of  the  algorithm  to  the  rotorcraft  de¬ 
scent  problem. 

Like  its  name,  the  algorithm  has  two  distinct  phases, 
the  gradient  phase  and  the  restoration  phase.  In  the  gra- 


dient  phase,  the  solution  is  modified  to  reduce  the  objective 
function.  In  the  restoration  phase,  the  constraint  error  is 
minimized  until  it  is  near  zero  within  a  specified  tolerance. 
For  example,  the  initial  conditions  for  certain  parameters 
may  be  simple  expressions  like  a  linear  variation  that  does 
not  satisfy  the  equations  of  motion.  The  restoration  phase 
adjusts  the  variables  until  all  of  the  constraints  are  satis¬ 
fied.  The  gradient  phase  also  introduces  constraint  error  as 
it  modifies  the  solution  to  minimize  the  objective  function. 
So  the  solution  alternates  between  gradient  steps,  which  im¬ 
prove  the  objective  function  but  introduce  constraint  errors, 
and  restoration  steps,  which  eliminate  the  constraint  error. 

To  ensure  a  stable  and  converging  solution,  a  limit  is 
placed  on  how  much  error  a  gradient  step  can  introduce.  If 
an  iteration  step  is  calculated  which  results  in  too  much  er¬ 
ror,  the  step  is  halved  until  the  resulting  error  is  within  lim¬ 
its.  Additionally,  the  objective  function  is  evaluated  after 
restoration  to  determine  if  the  gradient  step  was  successful. 
If  the  objective  function  is  larger  after  restoration  than  the 
solution  prior  to  the  gradient  step,  the  gradient  step  is  halved 
until  the  resulting  solution,  after  restoration,  is  improved. 

Initially,  the  differential  constraints  represented  equa¬ 
tions  of  motion,  the  non-differential  constraints  were  used 
to  limit  the  controls,  and  a  single  scalar  parameter  was  used 
to  scale  the  time  interval.  However,  additional  constraints 
were  required  and  increased  the  number  of  variables  in  the 
problem.  A  constraint  on  a  control  variable  adds  another 
control  variable.  A  constraint  on  a  state  variable  is  more 
costly,  adding  both  a  control  and  a  state  variable. 


Equations  of  Motion 


The  basic  equations  of  motion  for  horizontal  velocity,  verti¬ 
cal  velocity,  and  rotor  speed,  and  the  equations  for  rotor  in¬ 
flow  are  adapted  from  reference  1.  The  equations  of  motion 
are  obtained  from  the  sum  of  forces  in  the  x  and  z.  directions 
and  sum  of  torque  about  the  rotor  shaft.  After  modification 
to  change  the  analysis  variables  from  Cf  .  and  Cj,  to  Cj  and 
a  as  control  variables,  they  are  given  in  non-dimensional 
form  by 


d  X  z 

d 

ZTvx 

d  x 
d 

-j-® 

dx 


dx 


+ 


z  =  V, 


r®2  j  Qo 
a  I  8 
/  Cj/ o 


1+6 


CrY 

vj 


(1) 

(2) 


[1+4.6m2]  +  ^A|(3) 
(4) 


The  vertical  force  equilibrium  in  equation  1  is  the  sum 
of  the  weight  (the  equation  is  non-dimensionalized  and  nor¬ 
malized  by  Ct0 ,  so  the  1  represents  the  weight),  the  vertical 
component  of  thrust,  and  the  vertical  drag  of  the  fuselage. 
Note  that  the  weight  is  positive  while  the  thrust  and  drag 
terms  are  negative  because  the  sign  of  z  is  positive  down. 
Horizontal  equilibrium  in  equation  2  is  similar  without  the 
1  representing  non-dimensionalized  weight. 

The  torque  equilibrium  in  equation  3  is  a  sum  of  a  semi- 
empirical  expression  for  profile  power  and  the  momentum 
theory  approximation  of  induced  power.  Of  particular  note 
is  the  third  term  in  brackets.  This  term  is  intended  to  en¬ 
courage  the  optimization  algorithm  to  stay  below  stall  by 
rapidly  increasing  torque  above  a  specified  stall  thrust  coef¬ 
ficient.  A  typical  value  for  ns  is  20  or  30,  so  when  the  term 
in  parentheses  exceeds  unity,  the  torque  rises  very  quickly. 
Equation  4  is  only  needed  to  track  altitude  to  determine 
when  the  helicopter  reaches  the  ground. 

For  these  equations,  Cj  and  a  are  control  variables,  vz, 
vx,  (0,  and  z  are  state  variables,  and  the  other  variables  are 
either  specified  or  calculated  from  the  states  and  controls. 
X  is  calculated  with  corrections  for  vortex  ring  state  and 
ground  effect  based  on  the  horizontal  and  vertical  speeds 
and  the  non-dimensional  height  above  the  ground. 

The  rotor  induced  velocity  A  is  assumed  to  be  quasi¬ 
steady  and  is  given  by 

A  =  1.15Ahover///G  (5) 

This  equation  is  based  on  the  hover  induced  velocity 
Ahover  =  ^ / (2pAQ.^R2) .  The  1.15  factor  is  an  analysis  vari¬ 
able  but  is  not  normally  modified.  //  and  fa  are  factors 
to  correct  for  vortex  ring  state  and  ground  effect.  Like  the 
equations  of  motion,  they  were  taken  from  reference  1 ;  they 
originally  are  found  in  references  5  and  6,  respectively.  Far 
from  the  ground  and  at  normal  airspeeds,  both  factors  ap¬ 
proach  unity. 

The  inflow  curve  //  has  two  expressions  depending  on 
the  rotor  state.  A  transcendental,  but  analytic  momentum 
theory  expression  is  used  for  normal  forward  flight.  For 
nearly  vertical  flight,  where  momentum  theory  is  not  valid, 
an  empirical  expression  for  the  vortex  ring  state  is  used.  The 
equations  are  as  follows: 


//CM)  = 

/l/fV^+C x  +  fi)2)  ,  (2r+3)2+z2  >  1 

\x(0.373x2+0.598z2  —  1.991)  ,(2Jc+3)2+z2<  1 


The  parameters  x  and  z  are  normalized  vertical  and  hor¬ 
izontal  velocities, 


'fclimb  vx  sin  a  —  v-  cos  a 


Miover 


^hover^O^ 


(7) 


z  = 


'Mioriz 


vx  cos  a  +  v7  sin  a 


^hover  Ai;o\  e:  ““0  ^ 

Ground  effect  is  modeled  with  the  expression 


(8) 


,  =  1 _ (-vz  +  Acosa)2 _ 1_ 

G  (— v.  +  Acosa)2  +  (v'x  +  Asina)2  (4z)2 

where  z.  is  the  non-dimensional  rotor  height  above  the 
ground.  Note  that  the  normal  analysis  variable  z  is  the  skid 
or  wheel  height,  so  z  in  the  above  expression  must  be  cor¬ 
rected  by  the  height  of  the  rotor  with  the  helicopter  on  the 
ground.  Because  ground  effect  is  only  significant  near  the 
ground,  where  z  is  small,  the  distinction  is  important. 


Non-Dimensionalization  and  Scaling 

The  analysis  operates  on  the  equations  in  non-dimensional- 
ized  and  scaled  form.  The  non-dimensionalization  follows 
standard  convention  for  rotorcraft  for  the  state  variables. 
Additional  scale  factors  are  required  for  the  optimal  con¬ 
trol  procedure.  The  optimal  control  problem  is  well-scaled 
if  similar  changes  in  the  variables  produce  similar  changes 
in  the  objective  function.  For  example,  Cj  is  on  the  order  of 
0.005  whereas  forward  velocity,  even  non-dimensionalized 
by  rotor  speed  may  be  0.4.  Scale  factors  are  applied  to  bring 
the  order  of  magnitude  of  each  variable  to  approximately  1 . 
The  states,  controls,  and  scalar  parameter  appear  in  the  con¬ 
straint  equations,  therefore,  as  X/X,  U /U,  IT/fl.  Though 
these  factors  are  present  in  the  code,  for  clarity,  in  this  pa¬ 
per,  equations  are  presented  in  non-dimensional  form  with 
this  scaling  already  applied.  The  numerical  values  for  the 
scale  factors  are  given  later  when  all  of  the  states  and  con¬ 
trols  have  been  defined. 

The  terminal  time  is  T y,  as  yet  unspecified.  The  optimal 
control  algorithm  operates  on  a  specified  interval  from  0  to 
1.  So  an  additional  normalization  of  the  scaled,  dimension¬ 
less  time  is  required, 

%  =  f  do) 

Tf 

This  ensures  that  the  time  variable  the  analysis  uses 
ranges  from  0  to  1.  Because  non-dimensionalization,  a 
scale  factor,  and  the  normalization  have  been  applied,  dif¬ 
ferentiation  and  simply  getting  time  in  seconds  from  the 
unit  time  vector  in  the  analysis  can  be  confusing.  The  dif¬ 
ferential  operators  are  given  in  terms  of  f,  T,  and  the  single 
scalar  scale  factor  fli  by 


()V 


d  d 
d ^  Xf  dx 


V_  d 

floFli  dt 


(11) 


Table  1.  State  and  control  variables  used  in  the  baseline 
analysis 


Variable 

Symbol 

Description 

2ft 

vz 

Vertical  velocity 

x2 

Vx 

Horizontal  velocity 

*3 

(0 

Rotor  speed 

X4 

z 

Altitude 

Ul 

Cj 

Rotor  thrust  coefficient 

U2 

a 

Disk  plane  angle 

the  other  scale  factors  in  separate  program  units.  Each  of 
equations  1-4  therefore  contain  an  additional  if  factor  in 
the  code  to  reflect  the  ( )v  differentiation.  To  obtain  t  in 
seconds  from  the  unit  interval  time  c, , 


t  =  S 


1 

a0ni 


(12) 


For  the  remainder  of  the  paper,  the  if  will  be  explicitly 
present  as  the  additional  differential  equations  are  devel¬ 
oped. 

The  analysis  variables  in  terms  of  the  state  and  control 
vectors  X  and  U  are  given  in  Table  1 .  Additional  state  and 
control  variables  added  as  a  result  of  the  modifications  to 
the  baseline  model  are  discussed  below. 


Implementation  for  Power  Loss  Application 

The  preceding  discussion  describes  the  math  model  and  ba¬ 
sic  equations  of  motion  for  the  helicopter  in  flight  and  rep¬ 
resent  only  a  slight  departure  from  previous  work.  The  re¬ 
mainder  of  the  paper  is  the  implementation  and  application 
of  this  analysis  to  the  problem  of  descent  or  flyaway  with 
power  loss.  The  additional  constraints  for  the  problem  and 
boundary  values  are  discussed. 

At  several  points  in  the  course  of  development,  it  was 
observed  that  the  optimal  control  algorithm  would  produce 
flight  paths  that  violate  physics  but  are  well  within  the  con¬ 
straints  of  the  analysis  at  that  point.  It  is  at  the  same  time 
fascinating  and  frustrating  to  observe  the  innovative  solu¬ 
tions  that  are  produced  when  the  analysis  is  unleashed  on 
a  wide  range  of  initial  conditions.  It  is  fascinating  because 
some  solutions  are  not  necessarily  intuitive  to  a  human  an¬ 
alyst,  yet  are  perfectly  logical  when  the  actual  constraints 
or  lack  of  constraints  on  the  flight  path  are  considered.  But 
because  the  algorithm  will  only  obey  laws  of  physics  when 
forced,  each  innovative  yet  impossible  solution  mandates  a 
different  approach  or  additional  constraints  on  the  problem 
and  the  attendant  analysis  and  run  time  penalty  associated 
with  them. 


In  this  paper,  the  equations  are  given  in  if(d/di)  form 
since  it  is  the  easiest  to  understand.  The  computer  imple¬ 
mentation  specifies  equations  in  this  form  also,  and  applies 


Four  modifications  were  made  to  the  basic  model  and 
affect  both  controls  and  states.  Initially,  the  analysis  was 
for  power  off  descent  to  landing  only.  Three  changes  were 


implemented  at  that  time:  an  altitude  constraint,  rate  con¬ 
trols,  and  control  limits.  The  last  modification  was  for  par¬ 
tial  power,  which  was  added  in  conjunction  to  the  exit  and 
end  conditions  discussed  above.  These  four  modifications 
are  discussed  below. 

Altitude  Constraint 

The  first  modification  to  the  analysis  is  the  most  amusing 
and  serves  as  a  good  example  to  describe  how  state  con¬ 
straints  are  implemented  in  the  analysis.  The  source  of  the 
problem  was  the  observation  that  at  high  forward  speed  but 
low  altitude  initial  conditions,  the  optimal  flight  path  is  for 
the  helicopter  to  descend  below  ground  level  and  ascend  to 
the  landing  point  to  quickly  decelerate  and  meet  the  zero 
altitude  end  condition.  Because  the  analysis  was  only  con¬ 
strained  to  be  zero  at  the  final  time,  it  was  free  to  take  on 
negative  altitude  elsewhere  in  the  flight  path.  This  behavior 
is  replicated  in  nature  when  a  bird  approaches  the  edge  of 
a  building  or  a  wire.  It  approaches  slightly  below  the  alti¬ 
tude  of  the  landing  spot  and  as  it  nears  its  target,  ascends 
quickly  to  alight  on  the  roof  or  wire  with  small  horizontal 
and  vertical  speed.  While  it  is  encouraging  that  the  flight 
path  calculated  by  the  algorithm  is  found  in  nature,  it  is  not 
acceptable  for  the  purposes  of  autorotative  descent  of  a  he¬ 
licopter. 

A  constraint  was  added  to  force  the  analysis  to  keep 
the  helicopter  above  ground  prior  to  touchdown.  The  Se¬ 
quential  Gradient-Restoration  algorithm  operates  on  equal¬ 
ity  constraints,  so  some  additional  conversion  is  necessary 
to  implement  inequality  constraints,  such  as  (z  >  0).  Each 
inequality  constraint  must  be  converted  to  an  equality  con¬ 
straint  in  order  to  be  implemented  in  the  analysis.  The  pro¬ 
cess  adds  a  control  variable  and  an  equality  constraint  to  the 
analysis. 

Initially,  equation  4  tracks  the  vehicle  altitude.  Cor¬ 

rected  for  normalized  time,  the  equation  becomes 

;v  =  Tfvz  (13) 

z  is  positive  down,  so  it  ranges  from  (0  <  z  <  z/).  So  zo  =  0 
and  Zf  is  equal  to  the  initial  altitude.  The  desired  inequality 
constraint  therefore  is 

z<z/  (14) 

This  is  converted  to  an  equality  constraint  by  using  a  slack 
variable  z: 

z  — z/  +  z2  =  0  (15) 


zv  +  2zzv  =  0  (16) 

zv  =  xfUz=xfU5  (17) 

Now  equation  13  is  replaced  with  equation  17  in  the 
equations  of  motion  and  the  altitude  is  tracked  through  the 
slack  variable  z  rather  than  z  itself.  Though  generally  not 
so,  in  this  case  because  of  the  sign  conventions,  it  is  more 
convenient  to  track  z  since  it  is  zero  at  'Ey,  although  it  is 
initially  x/zo  rather  than  zo-  Altitude  is  easily  recovered 
through  equation  15.  After  substituting  equations  13  and 
17  into  equation  16,  a  non-differential  equality  constraint  is 
is  added  to  the  S  matrix  as 

vz  +  2zUz  =  0  (18) 

A  numerical  problem  arose  as  a  result  of  equation  18  in 
the  S  matrix.  Because  the  end  constraint  forces  z  to  zero  at 
the  final  time,  at  that  time  S  becomes  insensitive  to  Uz  at  that 
point.  At  times  UZf  would  diverge  to  large  values,  causing 
the  entire  problem  to  be  ill-conditioned.  A  static  offset  to  z 
was  applied  so  that  the  helicopter  is  instead  constrained  to 
remain  above  a  slightly  negative  altitude.  To  an  engineer¬ 
ing  approximation,  this  accomplishes  the  same  purpose  but 
avoids  the  numerical  implementation  issue. 

The  lower  altitude  bound  need  not  be  constant,  it  is  only 
required  to  be  a  continuously  differentiable  function.  A 
more  rigorous  method  is  envisioned  where  the  bound  on  the 
inequality  constraint  is  not  constant  but  a  line  which  crosses 
the  x-axis  between  the  T„,-i  and  x /  time  steps.  Such  a  con¬ 
straint  would  ensure  that  the  helicopter  remains  above  the 
surface  of  the  ground  at  all  times  except  'Ey  and  at  'Ey,  IJZ  is  a 
slightly  negative  number  so  a  proper  derivative  for  S%  can  be 
calculated.  Such  an  approach  has  not  yet  been  implemented 
or  tested. 

Rate  Controls 

The  approach  of  varying  Cj  and  a  as  controls  works  well 
when  the  initial  conditions  for  the  optimal  control  are  such 
that  there  is  plenty  of  altitude  and  forward  speed  to  per¬ 
form  the  autorotative  maneuver.  However,  for  cases  on 
the  fringe  and  in  regions  of  the  altitude/velocity  spectrum 
where  a  landing  is  not  possible,  many  iterations  are  per¬ 
formed  as  the  algorithm  tries  reduce  the  objective  function 
as  much  as  possible.  In  this  case,  often  small  variations 
in  the  controls  are  magnified.  The  Cj  and  a  controls  are 
discretized,  merely  a  list  of  values  at  each  of  the  nt  steps. 
There  is  nothing  preventing  large  step  changes  in  thrust  or 
disk  angle  from  one  time  step  to  the  next,  and  these  often 
occur  as  the  landing  task  becomes  more  difficult. 

Though  thrust  can  vary  relatively  quickly,  the  pitch  rate 
of  the  aircraft  has  a  definite  limit,  so  step  changes  in  a  are 
not  realistic.  Rate  controls  solve  this  problem.  Rather  than 


Taking  the  derivative  with  respect  to  tj , 


controlling  Cj  and  a  directly,  C'T  and  a*  are  controlled  and 
Cj  and  a  become  state  variables.  Step  changes  in  C*r  or  a* 
are  smoothed  by  the  integration  operation.  This  conversion 
replaces  the  two  controls  and  adds  two  additional  states  to 
the  analysis.  At  the  same  time,  each  rate  was  limited  by  an 
inequality  constraint,  which  added  two  additional  controls 
and  two  additional  non-differential  constraints  (one  control 
and  one  constraint  for  each  rate  limit). 

The  thrust  and  pitch  rates  should  vary  from  zero  to  a 
maximum  rate  (positive  or  negative).  This  can  be  written 
mathematically  as 


o  <(C*T)2<(C*Tmax)2  (19) 

0  <  (a*)2  <  «J2  (20) 


prevent  the  analysis  from  using  negative  thrust.  Some  care¬ 
ful  crafting  of  the  constraint  equation  provided  both  upper 
and  lower  limits  on  thrust  with  a  single  equation.  The  user 
inputs  are  the  maximum  and  minimum  thrust  coefficients 
allowed  for  the  helicopter.  The  maximum  is  unrelated  to 
the  stall  thrust  coefficient.  It  would  normally  be  higher  but 
can  be  the  same  or  lower  if  that  is  desired  by  the  user.  An 
additional  state  X-]  and  control  U(  are  added  as  a  result  of 
the  limit.  The  equality  constraint  is  given  by 


+  Cj 


(25) 


Taking  the  derivative,  the  equality  constraint  is  obtained 


The  inequality  constraints  were  converted  to  equality 
constraints  in  the  same  manner  as  described  previously.  The 
controls  were  implemented  as  U\  =  Cj  and  U2  =  (X*  being 
the  controls  for  thrust  rate  and  pitch  rate,  and  f/3  and  U4  as 
the  slack  variables  for  the  thrust  and  pitch  rate  limits,  re¬ 
spectively.  The  resulting  equations  for  the  S  matrix  are 

f/2+t/32-(c;max)2  =  0  (2i) 

Ui  +  ui-(a*mj2  =  0  (22) 


^ CT  -  Cr™n^Crmax^  Cj  +  CtCtV  =  0  (26) 

C?  =  yuCr  =  T fU6  (27) 

and  the  expression  for  Cj  is  given  in  equation  23.  Equation 
26  becomes  an  additional  non-differential  constraint  in  the 
S  matrix  and  equation  27  is  an  additional  differential  con¬ 
straint  in  the  <p  matrix. 

Partial  Power 


Additional  state  variables  A5  =  Ct  and  A),  =  a  were 
added,  resulting  in  additional  differential  constraints: 


xj  =  Cj  =  x fC*T  =  ZfUl  (23) 

Xj  =  av  =  x fa*  =  xfU2  (24) 

Even  if  the  optimal  control  procedure  prescribes  a  step 
change  in  Cj  or  a*,  the  time  histories  for  Cj  and  a  are 
continuous. 


Thrust  Limiting 

Although  the  torque  equation  (3)  provides  a  penalty  of 
rapidly  increasing  torque  above  a  specified  thrust,  explicit 
constraints  on  thrust  were  still  found  to  be  necessary.  Given 
unlimited  thrust  potential,  the  algorithm  would  produce 
control  schedules  which  resulted  in  a  Cj/o  greater  than  0.2 
and  sometimes  even  0.3  during  a  flare.  The  massive  torque 
increase  would  rapidly  slow  the  rotor  below  0.5f2o  instanta¬ 
neously  rather  than  cause  the  optimal  control  to  avoid  such 
high  thrust. 

An  explicit  thrust  limit  was  implemented  in  a  similar 
manner  to  that  for  the  altitude  limit.  In  this  case,  though,  in¬ 
dependent  upper  and  lower  limits  on  thrust  were  desired  to 


For  a  complete  power  loss,  the  torque  expression  in  equa¬ 
tion  3  is  set  equal  to  zero.  For  partial  power,  the  engine  can 
provide  torque  to  the  rotor.  For  a  real  engine,  a  control  sys¬ 
tem  or  person  adjusts  the  throttle  to  maintain  rotor  speed. 
If  the  torque  available  is  not  sufficient,  the  rotor  speed  will 
droop,  but  if  excess  torque  is  available,  the  control  system 
will  reduce  fuel  flow  to  maintain  100%  RPM.  This  level 
of  detail  is  not  replicated  in  the  current  analysis.  Rather,  a 
maximum  torque  is  known  from  user  input.  If  the  rotor  de¬ 
mands  more  than  the  maximum  available,  rotor  speed  will 
droop  like  a  real  aircraft,  but  if  excess  torque  is  available, 
the  algorithm  can  overspeed  the  rotor.  A  control  variable  Qe 
ranging  from  0  to  the  available  torque  specifies  how  much 
torque  is  applied  to  the  rotor  system.  A  second  (slack)  con¬ 
trol  variable  Uqb  is  also  required  to  convert  the  inequality 
constraint  to  an  equality  constraint. 

For  the  partial  power  case,  first  equation  3  is  modified  to 


d(0 

dx 


yco2 

a 


~Qe  + 


Cd0 
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1  + 
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(  Ct / 
V(Cr/a)i 


[1+4.6M2]  +  ^aJ(28) 


The  inequality  constraint  0  <  Qe<  go  is  converted  to  an 
equality  constraint  for  the  S  matrix  given  by 


Table  2.  State  and  control  variables  used  in  the  complete 
analysis 


Variable 

Symbol 

Description 

A, 

% 

Vertical  velocity  state 

a2 

Vx 

Horizontal  velocity  state 

a3 

CO 

Rotor  speed  state 

a4 

z. 

Slack  variable  (state)  for  altitude 

As 

Ct 

Rotor  thrust  coefficient  (state) 

a6 

a 

Disk  plane  angle  state 

a7 

CT 

Slack  variable  (state)  for  Cj 

t/i 

c* 

Thrust  rate  control 

u2 

a* 

Pitch  rate  control 

U3 

UC*T 

Slack  control  for  thrust  rate  limit 

Ua 

Ua* 

Slack  control  for  pitch  rate  limit 

U5 

Uz 

Slack  control  for  altitude 

u6 

UcT 

Slack  control  for  thrust  limit 

Vn 

Qe 

Engine  torque  control" 

Us 

Uq, 

Slack  control  for  torque  limit3 

flOnly  present  for  partial  power  analysis 


Q2e  +  Ul-Ql  =  0  (29) 

The  complete  set  of  states  and  controls  is  detailed  in  Ta¬ 
ble  2.  When  compared  to  Table  1,  it  is  clear  that  to  enforce 
limits  on  the  problem,  the  size  of  the  problem  grows  signifi¬ 
cantly.  In  its  current  state,  hard  limits  exist  on  every  control 
in  the  system  and  on  some  states.  The  original  states  vz,  vA- 
and  co  remain  unconstrained.  The  work  in  references  2  and 
3  imposed  a  limit  on  vertical  descent  rate,  but  such  a  limit 
has  not  yet  been  found  to  be  warranted.  The  nature  of  calcu¬ 
lating  the  optimal  path,  possibly  in  combination  with  thrust 
and  thrust  rate  limits,  seems  to  prevent  the  analysis  from  al¬ 
lowing  a  high  vertical  sink  rate.  Verification  and  validation 
efforts  in  the  future  may  reveal  that  a  hard  constraint  on  vz 
is  necessary,  but  results  to  date  suggest  it  is  more  likely  that 
co  will  require  limiting. 


End  Conditions  and  Initial  Path 

The  optimal  control  operates  on  an  initial  path  from  the  ini¬ 
tial  conditions  (altitude,  speed,  and  power  available)  to  end 
conditions  which  depend  on  the  power  available.  If  suffi¬ 
cient  power  is  available  to  fly  away,  the  analysis  attempts  to 
calculate  a  flyaway  path.  If  insufficient  power  to  fly  away  is 
available  or  the  flyaway  path  calculation  fails,  the  a  descent 
to  landing  is  calculated.  For  each  situation,  an  initial  path 
must  be  supplied  to  the  optimal  control  algorithm,  prefer¬ 
ably  one  that  is  close  to  the  final  path  and  satisfies  the  con¬ 
straints.  Obtaining  an  initial  path  close  to  the  final  path  is 
difficult,  after  all,  finding  the  final  path  is  the  purpose  of  the 
algorithm  in  the  first  place.  But  finding  a  path  that  meets  or 
nearly  meets  the  constraints  is  achievable. 

First  it  is  important  to  understand  the  initial  conditions 
and  end  conditions  before  discussing  the  initial  path  be¬ 


tween  them.  Only  steady,  level  flight  conditions  at  entry 
are  currently  used.  The  altitude,  airspeed,  and  power  avail¬ 
able  after  the  event  completely  define  the  initial  state.  The 
initial  thrust  coefficient  and  disk  angle  are  calculated  from 
the  gross  weight  and  initial  speed.  The  end  conditions  dif¬ 
fer  for  landing  and  mission  abort  flyaway,  but  are  straight¬ 
forward.  They  are  implemented  in  a  matrix,  y/  which  the 
optimal  control  algorithm  forces  to  zero  at  the  end  of  the 
time  interval.  Ideally,  the  initial  path  should  result  in  y/  =  0 
but  this  is  not  required  and  is  not  possible  in  some  cases  as 
discussed  below. 

For  landing,  the  aircraft  clearly  must  touch  the  ground, 
so  the  final  altitude  must  be  zero.  In  the  current  implemen¬ 
tation,  for  a  power-off  landing,  this  is  the  only  end  criterion. 
For  a  partial  power  landing,  an  additional  constraint  is  ap¬ 
plied  that  the  disk  angle  a  must  also  be  zero.  This  is  based 
on  an  assumption  that  with  partial  power,  the  landing  task 
is  easier  for  the  pilot  and  should  be  accomplished  without 
dragging  the  tail  on  the  ground. 

For  mission  abort  flyaway,  the  helicopter  must  be  in 
level  flight  (zero  sink  rate),  with  rotor  speed  at  100%,  and 
the  forward  velocity  must  be  the  minimum  power  speed. 
Mathematically  these  conditions  are  represented  by 


Vpower  off 


^partial  power 


V^flyaway 


m 


CO-  1 

vz 


Vx  -  Mopt 


(30) 

(31) 

(32) 


For  descent  to  landing,  the  initial  conditions  are  ob¬ 
tained  by  time  integrating  the  equations  of  motion  with  a 
complete  power  loss  and  no  change  in  controls.  Without 
power,  the  helicopter  cannot  sustain  steady  level  flight  and 
descends  to  impact  the  ground.  For  a  power  off  landing,  the 
resulting  flight  path  satisfies  the  constraint  error  to  numer¬ 
ical  precision  so  that  the  solution  can  proceed  immediately 
to  a  gradient  phase.  For  the  partial  power  landing,  the  con¬ 
straints  are  not  completely  met  and  a  restoration  is  required 
to  satisfy  them. 

Initial  conditions  for  flyaway  are  in  development.  As¬ 
signing  initial  conditions  in  this  case  is  not  straightforward 
because  there  is  not  a  discrete  event  such  as  striking  the 
ground  that  can  be  reached  by  time  integration.  For  fly¬ 
away,  the  pilot  must  transition  from  the  flight  condition  at 
the  power  loss  event  to  a  condition  where  the  helicopter  can 
sustain  level  flight  to  escape  without  striking  the  ground. 
Beyond  these  criteria,  there  is  no  requirement  to  judge  one 
flight  path  superior  to  another.  If  the  end  state  of  the  flight 
path  represents  a  steady  flyaway  condition,  the  analysis  has 
succeeded  in  showing  that  the  aircraft  can  fly  away.  Al¬ 
gorithmically,  if  the  restoration  satisfies  the  constraints  in 


equation  32,  there  is  no  need  for  gradient  steps  to  improve 
the  solution  and  the  analysis  exits. 

But  mathematically  some  initial  conditions  must  be  sup¬ 
plied  to  the  algorithm  to  begin  with,  particularly  the  interval 
time  T f.  Several  approaches  have  been  investigated  which 
work  well  in  some  instances  but  not  others.  The  first  is 
to  use  identical  initial  conditions  to  landing  and  simply  let 
the  restoration  phase  restore  to  a  flyaway  instead.  This  ap¬ 
proach  breaks  down  at  very  low  altitude  because  the  time 
interval  to  impact  for  a  complete  power  loss  is  very  short. 
If  the  flight  condition  is  at  very  high  or  low  speed  and  low 
altitude,  the  restoration  phase  cannot  transition  to  the  mini¬ 
mum  power  speed  in  such  a  short  time. 

A  second  option  is  to  make  all  states  constant,  as  though 
the  helicopter  continues  flying  with  sufficient  power  avail¬ 
able  for  a  specified  time.  When  the  power  is  removed,  the 
restoration  procedure  starts  with  a  flight  path  that  meets 
most  of  the  constraints.  Only  the  torque  and  vx  —  /ropt  con¬ 
straints  are  violated.  It  must  then  adjust  the  flight  path  to 
account  for  the  torque  loss  and  acceleration  or  deceleration. 
For  this  case,  there  is  still  the  issue  of  determining  a  suit¬ 
able  analysis  interval,  if  should  be  short  if  the  helicopter 
is  already  near  the  minimum  power  speed  but  longer  if  it  is 
far  above  or  below  it.  To  account  for  this,  the  time  required 
for  a  constant  acceleration  or  deceleration  to  the  minimum 
power  speed  is  calculated.  For  example,  if  the  initial  speed 
is  100  knots  and  the  minimum  power  speed  is  80  knots,  the 
time  for  a  constant  0.2g  deceleration  from  100  to  80  knots 
is  calculated  and  used  for  T /.  This  approach  also  works 
for  some  cases  but  not  others.  In  particular,  when  the  ini¬ 
tial  speed  high,  there  is  excess  energy  and  the  pilot  should 
be  able  to  gently  decelerate  in  level  flight  to  the  minimum 
power  speed  while  maintaining  altitude,  but  this  behavior  is 
not  always  replicated  by  the  analysis. 

Objective  Function  and  Exit  Criteria 

The  purpose  of  optimal  control  algorithm,  mathematically, 
is  to  minimize  an  objective  function.  The  objective  function 
consists  of  an  integral  part  and  a  non-integral  (scalar)  part, 

=  &  +  Sf.  The  most  important  physically  is  the  non¬ 
integral  part,  given  by 

&  =  Vl+wyxi  (33) 

This  represents  a  combination  of  the  final  vertical  and 
horizontal  components  of  velocity,  scaled  relative  to  each 
other  by  Wx.  The  scale  factor  determines  the  extent  the  op¬ 
timal  control  procedure  will  work  to  minimize  one  term  rel¬ 
ative  to  the  other.  The  terms  are  squared  so  that  they  are 
always  positive. 

The  integral  part  is  not  strictly  necessary,  but  has  been 
implemented  in  order  to  encourage  the  algorithm  rather 
than  formally  constraining  it.  The  integral  portion  of  the 
objective  function  is  currently 


&  =  lO(fi)  —  1)2(1  —  <^4) 

+  (V*-JU min)2(l -COs27I^)+Vj  (34) 

The  three  terms  are  obtained  from  engineering  judgment 
and  observation  of  solutions  the  optimal  control  algorithm 
produces.  The  first  term  penalizes  the  objective  function 
as  the  rotor  speed  deviates  from  100%  RPM.  The  (1  —  <f;4) 
causes  the  penalty  to  go  away  at  the  end  of  the  time  inter¬ 
val  where  the  rotor  speed  naturally  drops  during  the  flare. 
The  second  term  encourages  the  algorithm  to  approach  the 
minimum  power  speed.  It  also  has  a  time  factor  so  that  the 
term  diminishes  to  zero  at  the  beginning  and  end  of  the  time 
interval.  This  term  is  particularly  important  for  descents 
starting  from  hover.  The  mathematics  of  the  analysis  see  no 
difference  between  the  helicopter  going  backward  or  for¬ 
ward  from  a  hover  initial  condition.  With  this  term,  there  is 
a  mathematical  reason  to  go  forward.  The  final  term  penal¬ 
izes  the  vertical  velocity  at  all  times  in  the  time  interval  to 
provide  further  encouragement  to  the  analysis  to  not  build 
up  an  excessive  sink  rate. 

As  the  optimal  control  algorithm  minimizes  the  horizon¬ 
tal  and  vertical  speeds  at  impact,  a  solution  can  be  reached 
where  further  optimization  of  the  flight  path  will  not  change 
the  answer.  The  decision  between  an  attrition  and  a  forced 
landing  is  based  solely  on  vertical  impact  velocity  vCIit. 
Once  the  optimizer  reduces  the  vertical  impact  velocity  be¬ 
low  v'crit,  the  landing  is  classified  as  a  forced  landing  and 
further  optimization  is  not  necessary.  In  that  case,  if  the 
horizontal  velocity  criterion  vxi  <  W*Vcrit  has  been  met,  the 
optimal  control  procedure  exits  and  the  analysis  moves  on 
to  the  next  point.  If  the  horizontal  velocity  criterion  has  not 
been  met,  the  analysis  continues  until  it  is.  The  additional 
check  on  horizontal  velocity  is  required  to  account  for  high 
forward  speed  with  low  initial  altitude.  In  this  case,  the  ver¬ 
tical  sink  rate  may  be  less  than  vcr;t  after  the  initial  restora¬ 
tion  phase  while  the  horizontal  speed  is  very  high.  If  vx  is 
not  checked,  the  analysis  exits  at  this  point  and  incorrectly 
declares  the  point  a  forced  landing.  The  additional  vx  crite¬ 
rion  prevents  this  exit  and  allows  gradient  phases  to  occur 
to  minimize  the  horizontal  velocity  also. 

The  remaining  two  possible  criteria  normally  result  in  an 
attrition  because  the  optimal  control  algorithm  is  not  able  to 
reduce  the  impact  velocity  sufficiently  and  exits  for  one  of 
two  reasons.  First,  a  rolling  log  of  the  objective  function  af¬ 
ter  each  gradient  iteration  is  also  maintained.  If  the  objec¬ 
tive  function  has  not  been  reduced  by  1%  in  five  gradient 
iterations  (in  other  words,  for  iteration  i,  .4)  <  0.99J^_5), 
the  optimal  control  considers  the  current  solution  to  be  op¬ 
timal  and  exits.  Alternatively,  the  code  counts  gradient  iter¬ 
ations.  If  the  number  of  gradient  iterations  exceeds  a  preset 
amount,  the  solution  exits. 

There  is  a  third,  defunct  possibility  for  an  unsuccessful 
exit  from  the  optimal  control.  For  any  iterative  procedure. 


there  is  the  possibility  of  a  diverging  solution.  To  address 
that  situation  gracefully,  the  constraint  error  is  periodically 
checked  for  values  of  infinity  or  NaN  (not  a  number).  The 
constraint  error  is  calculated  from  all  of  the  variables  in  the 
analysis,  so  a  non-numeric  value  in  any  variable  naturally 
propagates  to  it.  If  such  a  value  is  detected,  a  flag  is  set 
which  causes  the  analysis  to  “fall  out”  of  any  nested  loops 
and  move  on  to  the  next  case. 

Validation  with  Autorotation  Data 

Autorotation  data  from  references  7  and  8  are  used  for  val¬ 
idation  of  the  power-off  performance.  The  data  are  from  a 
high  inertia  rotor  test  program  where  a  helicopter  was  out¬ 
fitted  with  additional  tip  weights  to  store  more  energy  in  the 
rotor.  The  effects  of  the  additional  energy  on  autorotation 
and  abrupt  maneuvers  was  then  evaluated. 

The  test  flights  involved  what  the  authors  called  “throt¬ 
tle  chops,”  where  the  engine  continues  to  run  but  is  rapidly 
reduced  to  idle  or  near  idle  speed.  A  temporary  provision 
was  added  to  the  DESCENT  analysis  to  simulate  this  ex¬ 
ponential  decay  in  throttle  by  allowing  the  torque  to  decay 
from  the  pre-event  level  to  a  near  zero  level  estimated  from 
the  test  data.  The  optimal  control  normally  uses  all  of  the 
torque  available  to  it  when  partial  power  is  present,  so  the 
torque  applied  to  the  rotor  in  the  simulation  is  similar  to  that 
in  the  test.  Although  achieving  a  successful  landing  in  the 
analysis  was  straightforward,  other  inputs  had  to  be  modi¬ 
fied  to  cause  the  flight  path  to  landing  to  closely  match  the 
test  data. 

To  obtain  the  appropriate  torque  parameters,  the  code 
was  modified  so  that  Qe  varies  over  the  time  interval. 
For  equation  29,  the  constant  Qq  is  replaced  by  a  time- 
dependent  ()(£,  ).  Initially,  it  is  equal  to  the  (calculated)  en¬ 
gine  torque  prior  to  the  power  loss.  Reference  8  provides 
torque  pressure  data  as  a  function  of  time.  It  is  not  reported 
how  torque  pressure  in  psi  relates  to  torque  in  ft-lb,  so  Qe 
was  scaled  so  that  that  the  ratio  of  initial  torque  to  the  final 
torque  was  the  same  as  the  ratio  of  the  initial  torque  pres¬ 
sure  to  the  idle  torque  pressure  from  the  test  data.  For  some 
cases  the  torque  pressure  at  idle  was  nearly  a  third  of  the 
pressure  prior  to  engine  loss,  and  for  others  it  was  nearly 
zero. 

The  first  case  to  be  considered  is  the  only  forward  flight 
case  from  reference  8.  It  is  for  the  baseline  rotor  flying  at  45 
knots  and  100  feet,  although  the  test  data  indicate  that  the 
initial  altitude  is  actually  114  feet.  For  this  flight,  the  test 
data  show  that  the  helicopter  lands  at  a  relatively  high  hori¬ 
zontal  speed  of  approximately  20-30  ft/sec  but  with  almost 
no  vertical  speed.  The  landings  were  made  to  grass,  and  the 
authors  describe  the  helicopter  sliding  to  a  stop  in  the  grass 
after  touchdown.  To  replicate  this  behavior,  the  weighting 
factor  Wx  was  set  to  0.05  and  Vm  to  1  ft/sec.  This  allows  a 
horizontal  speed  of  up  to  20  ft/sec  at  touchdown. 


Fig.  2.  Comparison  of  skid  height  from  DESCENT  anal¬ 
ysis  and  Reference  8  for  a  modified  OH-58A,  114  ft  and 
45  kt  initial  conditions. 

The  time  histories  of  skid  height,  horizontal  distance, 
and  rotor  speed  are  shown  in  Figures  2-4.  The  first  item 
to  verify  is  the  length  of  the  x-axis  in  the  plots.  As  de¬ 
scribed  earlier,  in  the  DESCENT  analysis,  time  is  non- 
dimensionalized  by  rotor  speed  and  then  normalized  so  that 
the  optimal  control  performs  its  operations  on  a  unit  time 
interval.  In  addition,  Tr  is  scaled  by  IT  i  in  the  computer 
program  though  those  scale  factors  have  been  left  out  of  the 
equations  in  this  paper  for  simplicity.  Considering  all  of 
that,  although  it  seems  trivial,  it  is  quite  important  to  verify 
that  the  time  scale  is  correct.  The  optimal  control  response 
time  is  not  exactly  the  same  as  the  test  data,  but  is  clearly 
indicating  seconds  and  not  revolutions  or  some  other  scaled 
quantity.  The  analysis  interval  came  out  almost  exactly  the 
same  with  a  Wx  setting  of  0.03,  but  the  horizontal  distance 
and  landing  speed  were  significantly  higher  in  that  case. 

The  altitude  time  histories  (Figure  2)  show  similar 
trends,  although  the  human  pilot  allows  altitude  to  decrease 
more  rapidly  and  then  floats  along  the  ground  at  about  10 
feet  for  a  few  seconds  before  touching  down.  The  optimal 
control  descends  gradually  and  touches  down  without  float¬ 
ing. 

The  plots  of  horizontal  distance,  shown  in  Figure  3,  are 
quite  similar.  The  touchdown  velocity  (as  indicated  by  the 
slope  of  the  curve  at  touchdown)  is  about  the  same,  though 
the  DESCENT  analysis  touches  down  earlier.  This  value 
as  well  as  the  total  distance  traveled  are  strongly  dependent 
on  Wx.  The  exact  values  can’t  be  “dialed  in”  with  Wx  be¬ 
cause  the  optimal  control  works  in  discrete  iterations  and  a 
single  additional  iteration  may  change  the  touchdown  speed 
significantly. 

Rotor  speed  is  shown  in  Figure  4.  Fike  the  horizon¬ 
tal  distance,  the  final  rotor  speed  is  consistent  with  the  test 
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Fig.  3.  Comparison  of  horizontal  distance  from  DE¬ 
SCENT  analysis  and  Reference  8  for  a  modified  OH- 
58A,  114  ft  and  45  kt  initial  conditions. 


Fig.  4.  Comparison  of  rotor  speed  from  DESCENT 
analysis  and  Reference  8  for  a  modified  OH-58A,  114 
ft  and  45  kt  initial  conditions. 

data,  despite  the  touchdown  time  being  earlier.  The  decay 
of  rotor  speed  is  very  sensitive  to  Crmax and  (Cj / G)s.  Be¬ 
cause  of  the  rapid  increase  in  torque  with  stall,  the  rotor 
speed  can  decay  very  quickly  if  the  rotor  is  allowed  to  gen¬ 
erate  thrust  significantly  above  stall.  The  analyst  does  not 
normally  have  the  luxury  of  test  data  to  adjust  these  param¬ 
eters,  but  estimates  of  them  should  be  known.  The  impor¬ 
tance  of  these  inputs  is  magnified  when  trying  to  match  a 
specific  flight  path  from  test  data.  For  normal  operation  of 
the  analysis,  a  range  of  flight  paths  may  be  possible  and 
only  the  impact  velocity  is  subsequently  used  for  the  Pk/d 


Fig.  5.  Comparison  of  skid  height  from  DESCENT  anal¬ 
ysis  and  Reference  8  for  a  modified  OH-58A,  27-ft  hover 
initial  condition. 

calculations. 

The  next  case  is  a  low  hover  case,  where  the  helicopter 
starts  from  a  hover  at  27  feet  at  the  time  of  the  engine  loss. 
For  this  case,  the  vertical  impact  velocity  in  the  test  data 
was  estimated  to  be  approximately  3  ft/sec,  so  vcl-it  was  set 
to  2.5  ft/sec.  If  set  to  3  ft/sec,  the  analysis  calculates  a  pure 
vertical  descent  and  does  not  enter  the  gradient  phase  at  all. 
A  vertical  descent  is  perfectly  reasonable,  but  a  pilot  would 
probably  initiate  some  forward  motion  if  for  no  other  reason 
than  to  see  where  he  was  going. 

However,  the  optimal  control  does  not  produce  an  an¬ 
swer  with  less  than  2.5  ft/sec  velocity  without  changing  pa¬ 
rameters.  The  reason  is  that  the  optimal  control  is  trying 
to  minimize  both  equations  33  and  34.  It  doesn’t  prioritize 
individual  parts  of  the  equations  unless  encouraged  to  do 
so  by  altering  the  relative  magnitude  of  different  terms.  In 
order  to  match  the  test  data,  equation  34  was  modified  so 
that  the  vx  term  was  divided  by  3.  This  encouraged  the  op¬ 
timal  control  to  place  comparatively  more  emphasis  on  the 
impact  velocities  than  forward  velocity.  Unmodified,  inter¬ 
mediate  solutions  with  impact  speeds  as  low  as  2.6  ft/sec 
are  obtained,  but  the  final  solution  after  the  program  exits 
with  50  gradient  iterations  has  an  impact  velocity  of  3.05 
ft/sec. 

The  results  from  both  the  modified  and  unmodified  re¬ 
sults  are  shown  in  Figures  5-7.  The  only  plot  where  there 
is  a  significant  difference  is  the  horizontal  distance.  Figure 
6.  The  DESCENT  solutions  are  virtually  identical  for  skid 
height  and  rotor  speed. 

The  skid  height  profile.  Figure  5,  matches  very  closely 
between  the  test  data  and  DESCENT  calculation.  Both  the 
modified  and  standard  DESCENT  calculations  are  shown. 


Fig.  6.  Comparison  of  horizontal  distance  from  DE¬ 
SCENT  analysis  and  Reference  8  for  a  modified  OH- 
58A,  27-ft  hover  initial  condition. 

The  calculated  flight  times  are  slightly  shorter,  about  3.8  sec 
vs  about  4. 1  sec  for  the  test  data.  Otherwise,  the  profiles  are 
very  similar.  The  slopes  of  the  curves  at  touchdown  are 
about  the  same,  indicating  that  the  3  ft/sec  approximation 
of  impact  velocity  in  the  test  data  was  accurate. 

The  horizontal  distance  profiles  are  shown  in  Figure  6. 
Here  there  are  significant  differences  between  the  modified 
and  standard  DESCENT  calculations.  One  difference  be¬ 
tween  the  calculations  and  test  data  is  that  the  test  pilot  ap¬ 
parently  made  no  attempt  to  decrease  horizontal  velocity  at 


Fig.  7.  Comparison  of  rotor  speed  from  DESCENT 
analysis  and  Reference  8  for  a  modified  OH-58A,  27-ft 
hover  initial  condition. 


touchdown.  The  optimal  control  reduces  horizontal  veloc¬ 
ity  to  1.6  ft/sec  in  the  modified  case,  and  1.0  ft/sec  in  the 
standard  case.  The  modified  case  travels  forward  approx¬ 
imately  9  feet  while  the  unmodified  calculation  is  nearly 
17.  These  bracket  the  test  data,  where  the  pilot  advanced 
12  feet.  This  illustrates  the  sensitivity  of  the  optimal  con¬ 
trol  to  the  objective  function  and  the  analysis  parameters. 
So  while  it  is  difficult  to  get  the  calculations  to  quantita¬ 
tively  match  the  test  data,  it  is  clear  from  the  time  histories 
of  skid  height  and  rotor  speed  that  the  code  is  producing  a 
qualitatively  similar  flight  path. 

Rotor  speed  is  shown  in  Figure  7.  Like  the  skid  height, 
the  rotor  speed  is  nearly  identical  for  the  two  calculations. 
The  calculated  rotor  speeds  at  touchdown  are  nearly  iden¬ 
tical  to  the  test  data,  with  only  the  slight  time  offset  be¬ 
tween  the  two.  The  profile  is  very  similar.  The  rotor  speed 
changes  very  little  for  a  second  or  two  and  then  the  pilot 
starts  the  flare,  bleeding  off  rotor  speed  much  more  quickly 
until  touchdown. 

The  final  case  is  a  high  hover  case,  where  the  helicopter 
starts  at  100  feet.  For  this  case,  shown  in  Figure  17  of  ref¬ 
erence  8,  the  mid-inertia  rotor  is  on  the  helicopter.  The  ad¬ 
ditional  inertia  causes  more  energy  to  be  stored  in  the  rotor. 
For  the  purposes  of  input,  the  only  change  to  the  helicopter 
model  was  changing  the  rotor  inertia  from  646  slug-ft2  to 
1100  slug-ft2.  For  this  comparison,  no  modification  to  the 
objective  function  was  made  in  the  DESCENT  analysis. 

The  results,  shown  in  Figures  8-10,  are  similar  to  those 
for  the  low  hover  case.  The  flight  time  is  about  one  second 
shorter,  7  seconds  vs  8  seconds  for  the  test  data.  The  calcu¬ 
lated  flight  path  has  significantly  more  horizontal  travel,  but 
the  rotor  speed  profile  is  similar.  Like  the  other  cases,  flight 
paths  with  horizontal  travel  from  as  little  as  40  feet  to  over 
200  feet  could  be  produced  with  small  changes  in  the  input. 

In  summary,  while  a  direct  numerical  comparison  is  very 
difficult,  a  qualitative  comparison  shows  that  the  analysis  is 
reproducing  the  trends  in  the  test  data  for  both  hover  and 
forward  flight  initial  conditions.  The  amount  of  energy  be¬ 
ing  dissipated  from  the  initial  potential  and  kinetic  energy 
of  the  helicopter  in  the  analysis  replicates  the  test  because 
the  velocities  at  the  end  of  the  flight  are  the  same,  and  the  ro¬ 
tor  energy  being  dissipated  is  the  same  as  shown  by  figures 
4,  7,  and  10.  For  the  purposes  of  calculating  Pk/d  values, 
it  is  sufficient  that  the  optimal  control  procedure  approxi¬ 
mates  the  energy  management  efforts  of  a  human  pilot. 

Partial  Power  Flyaway 

Although  still  in  the  development  stage,  the  partial  power 
flyaway  functionality  in  DESCENT  is  mature  enough  to 
produce  optimal  control  solutions.  An  example  is  presented 
here  to  illustrate  the  strengths  and  weaknesses  of  the  cur¬ 
rent  methodology.  For  this  case,  an  engine  loss  at  100  feet 
altitude  and  160  kts  forward  speed  is  simulated.  After  a 
single  restoration  phase,  a  flight  path  is  produced  where  the 


Fig.  8.  Comparison  of  skid  height  from  DESCENT 
analysis  and  Reference  8  for  a  modified  OH-58A,  100- 
ft  hover  initial  condition. 

aircraft  decelerates  to  the  minimum  power  speed  of  approx¬ 
imately  80  kts  in  about  17  seconds,  ending  at  100%  RPM 
and  no  sink  rate  as  specified  in  equation  32.  The  airspeed  is 
shown  in  Figure  1 1 . 

The  controls  are  shown  in  Figure  12.  Figure  12a  shows 
the  thmst  controls  while  Figure  12b  shows  the  pitch  con¬ 
trols.  In  these  plots,  the  rate  controls  are  shown  as  solid 
lines  and  the  resulting  displacements  are  shown  as  dashed 
lines.  The  rate  Cf  initially  increases  the  thrust  coefficient  to 
a  high  level  and  and  near  the  end  of  the  analysis  decreases  it 


Fig.  9.  Comparison  of  horizontal  distance  from  DE¬ 
SCENT  analysis  and  Reference  8  for  a  modified  OH- 
58A,  100-ft  hover  initial  condition. 


Fig.  10.  Comparison  of  rotor  speed  from  DESCENT 
analysis  and  Reference  8  for  a  modified  OH-58A,  100- 
ft  hover  initial  condition. 


Fig.  11.  Airspeed  for  partial  power  flyaway  transition 
to  minimum  power  speed. 

slightly.  For  this  case,  the  stall  thrust  coefficient  (Cj /o)s) 
is  0.14,  so  it  is  exceeded,  but  only  slightly.  Meanwhile,  the 
aircraft  pitches  back  to  about  20  deg  nose  up  to  achieve  the 
rapid  deceleration  in  Figure  11. 

The  effect  on  the  sink  rate  and  altitude  is  shown  in  Fig¬ 
ure  13.  Like  the  previous  plots,  the  rate  (v,)  is  shown  with 
a  solid  line  and  the  displacement  (altitude  £)  is  the  dashed 
line.  As  the  aircraft  immediately  pitches  back  and  increases 
thrust,  the  sink  rate  is  initially  negative.  As  the  helicopter 
bleeds  off  kinetic  energy  from  the  excess  airspeed,  it  climbs 
from  100  to  almost  140  feet  and  then  abruptly  drops  to  1 10 
feet  and  levels  out.  This  flight  path  may  seem  erratic,  but  it 
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Fig.  12.  Displacement  and  rate  controls  for  thrust  coeffi¬ 
cient  and  disk  angle  for  partial  power  flyaway  transition 
to  minimum  power  speed. 


is  important  to  note  that  no  optimization  has  occurred.  Only 
the  constraint  error  has  been  eliminated. 

The  effect  on  rotor  speed  is  shown  in  Figure  14.  The 
reason  for  the  high  thrust  coefficient  is  more  clear  consider¬ 
ing  that  the  rotor  speed  drops  to  about  65%.  A  significantly 
larger  thmst  coefficient  is  necessary  for  the  same  thrust  at 
such  a  low  rotor  speed.  But  the  decrease  in  thrust  coupled 
with  the  rapid  decrease  in  altitude  causes  the  rotor  speed  to 
recover  to  100%  at  the  end  of  the  analysis  interval  to  meet 
the  end  constraint.  In  reality,  a  pilot  would  never  allow  rotor 
speed  to  go  so  low,  and  probably  could  not  recover  it  if  it  did 
droop  that  much.  But  there  is  no  constraint  on  rotor  speed 
in  the  present  implementation,  so  such  large  variations  can 
occur  in  the  simulated  flight  path. 


Fig.  13.  Vertical  speed  and  altitude  for  partial  power 
flyaway  transition  to  minimum  power  speed. 


Fig.  14.  Non-dimensional  rotor  speed  for  partial  power 
flyaway  transition  to  minimum  power  speed. 


This  example  illustrates  how  the  optimal  control  meets 
all  of  the  constraints  placed  on  it,  and  how  if  it  is  not  suf¬ 
ficiently  constrained,  will  produce  a  flight  path  that  is  not 
physically  reasonable.  To  obtain  a  flight  path  more  repre¬ 
sentative  of  what  a  human  pilot  would  fly  will  require  mod¬ 
ifications  to  the  model.  There  are  two  paths  to  take.  Either 
additional  constraints  must  be  placed,  creating  more  burden 
for  the  single  restoration,  or  the  existing  constraints  can  be 
relaxed  and  their  effect  achieved  using  the  gradient  phase  to 
optimize  the  path.  In  the  latter  case,  the  gradient  phase  can 
improve  the  flight  path  over  many  gradient  steps  while  the 
current  implementation  only  allows  one  restoration  phase. 


Conclusions 

An  analysis  for  flyaway  or  descent  to  landing  with  power 
loss  has  been  developed  using  optimal  control  to  simulate 
the  pilot  reaction.  The  momentum  theory  vehicle  model  ex¬ 
tends  previous  work  by  using  rate  controls  and  additional 
constraints  on  the  optimal  control  procedure.  Based  on 
work  to  date,  the  following  conclusions  are  offered. 

1.  Descent  to  landing  and  partial  power  flyaway  have 
been  considered.  Suitable  initial  and  end  conditions 
for  landing  have  been  developed.  Partial  power  fly¬ 
away  is  still  in  development  and  will  require  addi¬ 
tional  refinement  to  consistently  obtain  satisfactory 
flight  paths. 

2.  For  benign  conditions,  simple  controls  such  as  hori¬ 
zontal  and  vertical  thrust  coefficients  or  thrust  coeffi¬ 
cient  and  disk  angle,  are  adequate  controls. 

3.  Rate  controls  and  additional  constraints  were  found 
to  adequately  address  specific  deficiencies  observed 
in  the  original  analysis  when  a  broad  range  of  initial 
speeds  and  altitudes  are  tested,  at  the  cost  of  additional 
analysis  complexity. 

4.  For  autorotation  from  a  forward  flight  initial  condition, 
the  optimal  control  procedure-determined  flight  path 
agrees  well  with  test  data  for  a  human-piloted  autoro¬ 
tation. 

5.  For  autorotation  from  hover,  the  analysis  replicates  the 
qualitative  trends  of  the  test  data.  By  adjusting  the  ob¬ 
jective  function  to  approximate  the  horizontal  touch¬ 
down  speed  in  the  test  data,  quantitative  agreement 
could  also  be  obtained. 
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